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ABSTRACT 

We consider diquark condensation in finite density lattice SU(2). We first 
present an extension of Vafa- Witten result, on spontaneous breaking of vector- 
like global symmetries, that allows us to formulate a no-go theorem for di- 
quark condensation in a region of the chemical potential-mass parameter 
space. We then describe a new technique to calculate diquark condensation 
at any number of flavours directly at zero external source without using any 
potentially dangerous extrapolation procedure. We apply it to the strong 
coupling limit and find compelling evidences for a second order phase transi- 
tion, where a diquark condensate appears, as well as quantitative agreement 
between lattice results and low-energy effective Lagrangian calculations. 
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1 Introduction 



Recently the expected scenario for the phase diagram of QCD in the chemical 
potential-temperature plane has changed: besides the hadronic and quark- 
gluon plasma phases, the existence of a new state of matter has been proposed 
by several groups [|I|. This new phase is characteristic of the high density, low 
temperature regime: the asymptotic freedom of QCD and the known insta- 
bility of large Fermi spheres, in presence of (whatever weak) attractive forces, 
results in a pairing of quarks with momenta near the Fermi surface (in anal- 
ogy with the Cooper pairing in solid state systems at low temperature). A 
condensation of quark pairs should be the distinctive signal of the new phase 
and has been indeed predicted using simplified phenomenological models of 
QCD. 

Unfortunately the lattice approach, the most powerful tool to perform 
first principles, non perturbative studies, is affected in the case of finite den- 
sity QCD by the well known sign problem that has prevented until now any 
step towards the understanding of this new phase. This is not the case for 
the SU(2) theory where the fermions are in the pseudo-real representation 
and finite density numerical simulations are feasible. 

In this paper we present a detailed study of diquark condensation for the 
unquenched two colours model. This work follows a previous one where we 
considered chiral and diquark susceptibilities and found strong evidence for 
a phase transition separating the ordinary low density phase from a high 
density one where the chiral condensate vanishes and the baryon number 
symmetry is broken^]. In the present work we mainly focus on an approach, 
based on the analysis of the probability distribution function (p.d.f.) of 
the order parameter for Grassmann fields 0, to compute the value of the 
order parameter directly at zero external source without using any potentially 
dangerous extrapolation procedure, mandatory in the standard approach. 

This technique has also another relevant feature. Since configurations are 
generated directly at zero diquark source, the sign ambiguity in the definition 
of the fermionic partition function is avoided. It turns out that we are not 
restricted to = 8 but any number of flavours can be in principle simulated. 

The study of diquark condensation in two colours QCD is an interesting 
topic by itself and we can hope to use some of the results to get insights 
about the three colours case. We are interested in the similarities between 
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the phase diagram of SU(3) and SU(2), however we have to be aware of the 
differences between the two models: 

• The quark-quark condensate {qq) is coloured for the N^. = ^ theory and 

colourless for = 2. This implies that in the SU(3) case the diquark 
condensation has to be interpreted as a Higgs-like mechanism leading 
to the breaking of the local colour symmetry, while in the SU(2) case 
we have spontaneous symmetry breaking (SSB) of the Ub{^) global 
symmetry associated with conservation of baryon number. 

• The zero temperature critical chemical potential, related to the mass of 

the lightest baryonic state, is different in the two cases: it is expected 
to be 1/3 of the nucleon mass for SU(3) and 1/2 of the pion mass for 
SU(2) (i.e. zero at vanishing quark mass). For Nc = 2, in order to have 
a phase diagram similar to that of SU(3), it is appropriate to consider 
a non zero bare quark mass. 

• The lightest baryonic state is a fermion for the tree colours model and a 

boson for the two colours one. 

Our results can be used to do a highly non trivial, quantitative check of 
the continuum predictions for the SU(2) theory, as for instance those from 
low energy effective Lagrangian calculations 0. 

In the next section we present an extension, to SU(2) gauge theory at 
finite chemical potential, of the Vafa-Witten theorem on the impossibility to 
break spontaneously baryon number conservation in a vector-like theory [^], 
finding an analytical bound, useful to check approximated models. In section 
III we give a review of the p.d.f. formalism for Grassmann degrees of freedom, 
emphasizing on the details of its application to the analysis of the order 
parameter for the U(l) symmetry associated to baryon number conservation. 
In section IV we present some details of the numerical algorithm. The last 
section is devoted to the presentation and discussion of numerical results. 

2 The Vafa-Witten theorem at finite fi 

In their original paper Vafa and Witten proved that the two point corre- 
lation function of any operator with non zero isospin or baryon number (or 
any other conserved quantity carried by fermions but not by gauge bosons) 
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falls off exponentially at large distances. From this result they concluded 
that any order parameter for such vector-like global symmetries is zero and 
that massless bound states can only be built by zero bare mass constituents. 

The two main ingredients that are necessary to the proof are (i) the anti- 
hermiticity of the massless Dirac operator (to bound the fermion propagator 
for a given gauge background) and (ii) the positivity of the integration mea- 
sure (used to extend the above bound to the averaged correlation function). 
These conditions limit the applicability of the theorem to vector-like theo- 
ries (Yukawa coupling generally invalidates the argument) at 6' = and zero 
chemical potential. 

We are here interested to the case in which the CP-violating vacuum an- 
gle vanishes but the chemical potential /i is non zero. While for SU(3) this 
breaks both conditions, for two colours QCD only condition (i) is not ful- 
filled. We will show how the theorem can be extended to the latter case for a 
range of values of fi that depends on the bare quark mass. Such new version 
of the theorem will allow us to exclude spontaneous breaking of baryon num- 
ber (i.e. the formation of a diquark condensate that breaks baryon number 
conservation) for some region of the parameter space. 

In this work we consider the theory regularized on the lattice. In such a 
case there is no need to consider smeared operators as in the original paper 
and we have to prove a uniform upper bound for the fermion propagator 

\{x\A~\fx)\0)\<aexpi-P\x\) (1) 

where A(/i) is the Dirac operator on the lattice in presence of a chemical 
potential term and a, /? are constants. A(/i) can be written as the sum 
of the Dirac operator at zero chemical potential (A(0) = iA + m) plus a 
term containing all the /i dependence: A(/i) = A(0) + e(/i). The matrix 
e(/i) = {e^ — 1)G + (1 — e^^)G'^ contains only the forward (G) and backward 
(G^) temporal links and is bounded: ||e(/i)|| < 2sinh/x. Expanding the 
inverse of the Dirac operator^ 

A-^(/x) = (/ + A-i(0)6(/i))-iA-i(0) 

= A-\0) - A-i(0)e(/i)A-i(0) + A-i(0)e(/i) A-i(0)e(/i) A-i(O) 

+ ••• (2) 

^The expansion is valid if || A^^(0)e(/x) || < 1. Using ||A^^(0)|| < and the bound 
on e(/i) the above condition is fulfilled if m > 2smhfi. At the end we will check that our 
result is inside this region. 
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and using an integral representation for A ^(0) we get the following series 
expansion for the fermion propagator: 

(x|A-i(/i)|0) = / rfroe-'""°(x|e-*"''^|0) 

JO 

p nOO 

+■■■ (3) 

The expectation values in the r.h.s. of (|) can be bounded depending on the 
value of the integration variables: since 

P„(x, 0) = Kx|e-«\(/i) . . . e-"^|0)| < ||e(M)r E ^(^o + ■ ■ ■ + r^)' 

k>x 

we get the relations 

^n[x,U) s I otherwise 
If we take the modulus of expression (|^) and use relations (H) we get 

\{x\A-\^)\0)\<J2qn + T.^n (5) 
n 

where 



X 



and Dn is the set of points that satisfy the condition ro+. . .Tn > x exp(— 2 
Clearly J2 Qn is bounded by a geometric series 

< e-ll^ll^' ^ if ||e(/i)|| < m (6) 

We have to find an (exponential) bound for the other term in the r.h.s. 
of (^). Since 

Iq= f rfroe"™"" = le-'"^"""""' 
JDo m 
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it is possible to calculate 



4= / t/ro...r„e-'"(^°+-^") 

JDr, 



using the change of variables 



To + Ti + . . . r„ = to 
Ti + . . . r„ = ti 



It is easy to check that 
and finally 



n\ dm"- 

^.„<|(_i,.m'£^ (T) 

Note that Jo is an analytic function of m for m 7^ so the r.h.s. of (0) 
is equal to the Taylor expansion of Io{m — ||e(/x)||) around Io{m) provided 
m > \\e{fi)\\. 

Putting together the last result with (|^) and the condition for the Taylor 
expansion in (0) we get the final bound 

\{x\A-\fi)\0)\ < ^_ ^1^^^^^^ (e-ll^ll" + e-(-"ll^('^)ll)--P(-2||A|l)) (g) 

valid for any m > ||e(yu)||. Since ||e(/i)|| < 2sinh/i we conclude that the 
fermionic correlation function decays exponentially at least for any 

m > 2 sinh fi (9) 

Note that in the fi limit we recover the original result of Vafa and Witten 
i.e. the propagator goes to zero exponentially for any nonzero mass. 

The formation of a diquark condensate, as follows from (^), is excluded for 
fi < jl = sinh~^(m/2). Near the chiral limit, jl is smaller than half the pion 
mass so condition (|^) is fulfilled for the expected critical point /ic = m.„/2. 
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Our result is interesting because it poses a rigorous limit to diquark con- 
densation and this limit can be used as a check to phenomenological ap- 
proximations to the model. It also shows how the Vafa-Witten theorem on 
spontaneous breaking of vector-like global symmetries can be extended, for 
some region of the parameter space, to a theory at non zero chemical poten- 
tial. 



3 The probability distribution function of the 
diquark condensate 



The use of the p.d.f. to analyze the spontaneous symmetry breaking in spin 
systems or Quantum Field Theories with bosonic degrees of freedom is a 
standard procedure. Less standard is its application to QFT with Grass- 
mann fields where the fermionic degrees of freedom have to be integrated 
analytically. 

The version of this method we use has been developed to extract the chiral 
condensate in the chiral limit from simulations of QFT with fermions|p, and 
will be used to study the vacuum structure of two colours QCD at non zero 
density and specifically to extract the diquark condensate at j = 0. We 
refer to the original paper for a full description of the p.d.f. technique 
and present a brief introduction focusing on the peculiarities of the diquark 
condensate case. 

In analogy with the study of chiral symmetry, we can construct a two 
component diquark condensate vector {{il)T2ip+ipT2ip).,i{ipT2ip—ipT2ip)) which 
transforms under U(1)b as a vector under rotations in the plane. Therefore 
we can take any of these two diquark condensates as order parameter c for the 
U(1)b symmetry associated to baryon number conservation. In the following 
we will use the first component. Then let a be an index which characterizes 
all possible (degenerate) vacuum states and Wa the probability to get the 
vacuum state a when choosing randomly an equilibrium state. If Cq, is the 
value taken by the order parameter in the a state, we can write 



where V is the lattice volume and the sum is over all lattice points. -P(c), 




X 



(10) 
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the p.d.f. of the diquark order parameter c, will be given by 



a 

= Mm ^ f [dU] [#] [d^]e-^oiu)+^A^^^}^ ^r2ij{x) + ^r2^(x) - c) 

V^oo Z J V ~^ 

The main point is that while -P(c) is not directly accessible with a nu- 
merical simulation its Fourier transform 

P{q) = Jdc e'^'Pic) (11) 

can be easily computed. Inserting in (|Tl|) the definition of -P(c) and using an 
integral representation for the 5-function, we can compute the integral over 
the Grassmann variables: 



^'4>A'4>+ig/V(4>T2'ip+'4>T2ip) _ pff^(^ 

where we have defined 



PfB(^) (12) 



-0) ^ ( ) . (-) 

and A is the usual lattice Dirac operator (it contains the mass and /x depen- 
dence). 

After some algebra we obtain: 

PviQ) = y [dU]e-'oiu)^MSl^d,tAfl (14) 

where Pvio) is the Fourier transformed p.d.f. of the diquark order param- 
eter at zero external source and finite volume, Z is the standard partition 
function. 



To take advantage of ([l^ we should be able to compute correctly the 
Pfaffian involved. This indeed turns out to be easy. The only ambiguity is 
related to the sign of PfB{y) = iy^det B{iq/V). In the next section we 
will analize some properties of the matrix B and show how it is possible to 
compute PfB{iq/V) once we have the eigenvalues of -B(O). 

As we will discuss in the next section, in standard simulations it is prac- 
tically impossible to determine the sign of the Pfaffians involved and one is 
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forced to choose Nj equal to a multiple of 8. A remarkable property of our 
approach is that, since we are able to compute PfB{iq/V), any value of Nf 
can be considered in the simulations. 

Once we have the Pvil) and hence, by simple Fourier transform, the 
p.d.f. of the order parameter we need to extract the correct value for the 
order parameter. To do that we first have to recall that diquark condensate is 
the order parameter of the Ub{^) symmetry associated with baryon number 
conservation and that {ipT2ip + ipT2ip) and i{il!T2ip — ipT2ip) are the components 
of a vector (in a plane) which rotates by an angle 2a when we do a global 
phase transformation of parameter a on the fermionic fields. Therefore, if 
Co is the vacuum expectation value of the diquark condensate (in the infinite 
volume limit) corresponding to the a-vacuum selected by a diquark source 
term after taking the zero source limit, -P(c) can be computed as 

P(c) = ^ f daS{c- Co cos(2a)) (15) 

which gives P(c) = l/(7r(cQ — c^)^/^) for — co < c < Co and P(c) = 
otherwise [^] . In the symmetric phase Co = and -P(c) reduces to a 5-function 
in the origin. 

The above results are valid in the thermodynamic limit while, at finite vol- 
ume, the non analyticities of the p.d.f. are absent. Even without performing 
a detailed finite size scaling analysis, we expect, for the finite volume p.d.f. 
Pv{c), a function peaked in the origin in the symmetric phase and peaked at 
some non zero value in the broken phase. This is indeed the behaviour we 
can observe in fig. ^ where the the p.d.f. of the smallest volume simulated 
(see Section 5) is reported at different values of fi. It is clear that, increas- 
ing the chemical potential, the vacuum starts to be degenerate signalling a 
spontaneous breaking of the baryon number conservation. 

We can also compare -Py(c) for two lattice volumes in the symmetric and 
broken phase. This is done in fig. 2 where we see clearly as, increasing the 
volume, the peak of the p.d.f. becomes sharper. To determine the value of 
the diquark condensate we used the position of the peak: a definition that 
clearly converges to the correct value in the thermodynamic limit. 

From fig. |^ we see also that data of the larger volume are more noisy 
(indeed we also get negative values for the p.d.f. in the broken phase). To 
have an estimate of the errors we calculated several distribution functions 
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for the 6^ lattice using the jacknife procedure. We saw that, except for the 
critical region, the position of the peak was very stable. 



4 Simulation scheme 

The standard way to study SSB is to introduce first an explicit symmetry 
breaking term in the action. If we do that for the diquark in the SU(2) model 
we have to add a term|^ j{'ipT2tp + ipT2ip) and, after integrating the Grass- 
mann field, the fermionic contribution for Nj = 4 quark flavours becomes 
proportional to the Pfaffian of a 4:V x AV matrix p|: 



where B{j) is deflned in the previous section and j is real. 

Using the relation T2AT2 = A* we can easily prove that -B(O) is antiher- 
mitian and det B[j) > for any j. It can also be shown that the eigenvalues 
of -B(O) are doubly degenerate and B^{0) is block diagonal with two hermi- 
tian blocks on the diagonal having the same eigenvalues. It follows that to 
compute det-B(j) for any j (i.e. to obtain all the eigenvalues of -B(O)) it is 
sufficient to diagonalize only one block of 5^(0) (reducing the problem to 
the diagonalization of a 2V x 2V hermitian matrix) and then take the two 
(imaginary) square roots of its (real and negative) eigenvalues. 

To avoid the sign ambiguity in (0) it is customary to consider a theory 
with Nf = 8 quark flavours where the fermionic partition function becomes 
ZfermU) = det-B(j). This hmitation can be overcomeby exploiting our abil- 
ity to work directly at zero diquark source. In the j = limit the sign ambi- 
guity disappears and the Pfaffian is positive definite since PfB{0) = det A 
and the last quantity is real and positive for any value of /i. Then we can 
easily consider any value of Nf writing Zferm{j = 0) = (det -8(0))^-^/^. 

If we are interested in the p.d.f. of the diquark order parameter at zero 
external source, the Pfaffian in ( pi^D can also be easily computed. The non de- 
generate eigenvalues of -B(O) can be written as ±iA„ {n = 1, ■ ■ ■ ,V) with real 
and positive A„. Using the relation ([T6|) we arrive at the following expression 




(16) 



^ n=l 




(17) 
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The sign ambiguity is solved noticing that ( |T7| ) has to be positive at j = 0. 
Increasing q, (|l^ changes sign each time q/V is equal to one of the A„, 
except possibly in case of degeneracy in the eigenvalues of B{0) resulting 
in the Pfaffian not crossing zero but tangent to the horizontal axis. This 
situation never occurred in our simulations. 

The procedure we presented in previous section can be used for any value 
of the gauge coupling but we have studied the phase structure of the theory 
in the limit of infinite gauge coupling (/? = 0). The main reason lies in 
the possibility to check our results with standard numerical simulations as 
well as analytical calculations. In this way we can restrict our efforts to the 
exploration of the phase space in the mass-chemical potential plane. 

To simulate the /3 = limit of the theory we have measured fermionic 
observables on gauge configurations generated randomly, i.e. with only the 
Haar measure of the gauge group as a weight. This choice implies a Gaussian 
distribution of the plaquette energy around zero which, according to the re- 
sults of Morrison and HandsQ, has a net overlap with the importance sample 
of gauge configurations at the values of fi and m used in our calculations. The 
validity of this procedure has also been tested comparing different physical 
observables (number density and chiral condensate) with Hybrid Montecarlo 
results^. 

We have considered the theory in a 4^ and 6"^ lattice diagonalizing 300 
gauge configurations in the first lattice volume and 100 in the second one. As 
we pointed out in the introduction the simulations have been performed at 
non zero quark mass, in order to study a physical situation closer to SU(3). 
We choosed m = 0.025, 0.05, 0.20 and values of the chemical potential ranging 
from fi = to fi = 1.0. The values of the diquark condensate are determined 
from the position of the maximum of the probability distribution function 
and the errors have been determined using a jacknife procedure. 

All numerical simulations have been performed on a cluster of PCs at the 
INFN Gran Sasso National Laboratory. 



5 Results and conclusions 

Here we present the results for the diquark condensate at j = as a function 
of the chemical potential. 

Fig. 3 contains a comparison of 6^ and 4^ results for Nf = 1 and m = 0.2. 
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We can clearly distinguish two symmetric phases separated by a broken one 
and two, possibly continuous, transition points. The high density symmetric 
phase has no physical relevance, being consequence of the saturation of all 
lattice sites with quarks. This phenomenum has nothing to do with contin- 
uum physics, being a pure lattice artifact. The physically interesting phase 
transition, i.e. the transition that has a continuum counterpart, is the first 
one 0] and it is the only one we will consider in the following. 

With only two available volumes we do not have the possibility to perform 
a serious finite size scaling analisys. Anyway it is evident that increasing 
the lattice volume the behaviour near the (physical) critical point is more 
singular. A similar qualitative result holds for all the available masses and 
flavour numbers thus strongly supporting the picture of singular behaviour 
in the infinite volume limit. 

We have also taken advantage of our simulation scheme considering dif- 
ferent values of flavours. Increasing Nf our operators become more noisy, 
especially near the transition point, but we can safely extend our calcula- 
tions up to Nf = 4. In fig. 4 we plot the diquark condensate for the largest 
mass and Nf = 1, 2, 4. We see clearly that the three data set are almost co- 
incident and we conclude that no dependence on Nf is evident. Once again 
this result is valid also for the smaller masses. 

The existence of spontaneous symmetry breaking with a non vanishing 
diquark condensate has already been predicted by Mean Field calculations [|| 
as well as numerical calculations H, ^ and effective models 0. We can use 
our relatively large data set to take a step forward and make a quantitative 
comparison with available analitical predictions. 

In previous work we used low energy effective Lagrangian calculations to 
check J 7^ results finding remarkable agreement 0. This motivated us to 
repeat the same procedure for our j = results. In this case the authors of 
[0] provide the following formula for the diquark condensate: 

[ otherwise 

where {iI!iIj)q, are the chiral condensate and the pion mass at zero chemical 
potential. 

We used ( |T8| ) to fit our data for the larger lattice near the critical point: 
plotting ^*{iljip)^ versus fi^ we expect a hnear dependence. We used all the 
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available masses and the Nj = 1 results that have smaller errors. In fig. 
5 the m = 0.025 case is considered showing both the data points and the 
resulting linear fit. In all cases the is good (see table 1) and we conclude 
that ([TSD gives a good description of the numerical results. 

In fig. 6 we have reported the diquark condensate as a function of /x for 
all the available masses. Clearly, away from the critical point, we can appre- 
ciate the dramatic effect of lattice discretization that prevents the diquark 
condensate to stay at the fi = chiral condensate value ([T8|). 

In Table 1 we compare the results of the fits with the /z = determinations 
of (ipip) and performed on a 6^ x 12 lattice at /? = 0. The determination 

Table 1: Parameters for the low energy effective Lagrangian predictions. 



m 


im^ (fit) 


2 


W)o (fit) 




d.o. f. 


0.025 


0.165(2) 


0.1696(11) 


1.325(10) 


1.31(2) 


1.1 


0.05 


0.236(5) 


0.2405(9) 


1.29(2) 


1.29(2) 


2.2 


0.2 


0.484(4) 


0.4841(7) 


1.438(13) 


1.23(2) 


3.3 



of the critical point and chiral condensate is in good agreement in the two 
data set, at least for the two smaller masses. Increasing m the becomes 
larger and the determination of the chiral condensate is less accurate. This 
however is not surprising since we expect the validity range of the low energy 
Lagrangian prediction (|18]) to be restricted to the small quark mass region 
where the gap between the pion mass and the first non-Goldstone excitation 
is large 0. 

What is more surprising is that, like in the j 7^ case [0, our P = 
calculations have a incredibly good agreement with a continuum (even if 
based on an effective model) prediction. Since the analytical predictions are, 
for the values of /i presented, well inside the validity region of the low energy 
approximation we can conclude that this gives indication of a very small 
dependence of lattice results on p. 

It would be very interesting to test this prediction by performing finite 
coupling calculations. This can be done using our approach to extract di- 
rectly the diquark order parameter at zero external source together with a 
HMC algorithm for generating configurations at /3 7^ 0. In this case we can 
compute Pv{(l) in (14) considering the ratio P f B{iq/V) / det A as an observ- 



13 



able which, for each gauge configuration and q < V, is a number of order 
one (or smaller). This program has to face two difficulties: at zero tempera- 
ture we need to increase our lattice temporal extent and correspondingly the 
computer time; on the other hand only unconclusive results are available for 
the II — thermodynamics of SU(2). 
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Figure 5: vs /i^ for 6^ lattice, m — 0.025 and linear best fit. 
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Figure 6: Diquark condensate for 6'^ lattice, Nf = 1, m = 0.025 (diamonds), 
0.05 (squares), 0.2 (stars) and corresponding fits using expression (piSl) . 
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